Configuration mixing description of the nucleus ^^S 
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We study the structure of the neutron rich '^^S isotope with modern configuration mixing methods 
based on the Gogny interaction including beyond- mean- field effects. Restoration of particle number 
and rotational symmetries are taken into account as well as shape mixing in the whole triaxial (/32, 7) 
plane. We obtain a qualitative agreement between the calculated spectrum and the experimental 
data reported recently. Very extended collective wave functions in the (/32,7) plane are found for 
the lowest states, corresponding to very flat potentials, indicating strong configuration mixing and 
supporting the weakening of the N = 28 magic number. 



Degradation and/or appearance of shell closures in nu- 
clei beyond the valley of stability is one of the main topics 
in nuclear structure for the last two decades or more. In 
particular, the magic number = 28 associated to the 
shell closure of a spherical harmonic oscillator plus spin- 
orbit single particle potential has been extensively stud- 
ied in neutron rich nuclei both theoretical and experimen- 
tally. For instance, the low excitation energy measured 
for the first 2+ state in the ^^Si isotope has determined 
the removal of such a shell closure in this case [1]. Indeed, 
the nucleus ^^S is in between the spherical doubly-magic 
^^Ca and the oblate deformed ^^Si isotopes. Therefore, 
coexistence of different shapes and/or configurations are 
expected to occur in this nucleus. Its transitional char- 
acter has been experimentally established by the appear- 
ance of a low energy O2" excited state [2] and the small 
value of the E{2f) and the large reduced transition prob- 
ability B{E2,2f of) [3]. Very recently, the values of 
the monopole strength p'^{E0^02 Of) and the transi- 
tion B{E2,2f O2) have been measured [4]. Based on 
shell model calculations and a simplistic two-level mixing 
model, these new data have been interpreted as an indica- 
tion of a prolate (neutron 2p-2h) -spherical (Op-Oh) shape 
coexistence. Additionally, some new levels have been un- 
ambiguously identified in Ref. [5], specially the 4^ level 
which mainly decays into the 2^ state. According to 
shell model calculations, this state would correspond to 
the band head for a K = 4 rotational band (neutron Ip- 
Ih configuration) although the 43" state with the neutron 
2p-2h configuration of the prolate band appears almost 
degenerated in the theoretical predictions. 
On the other hand, both mean field calculations with 
Skyrme and Relativistic interactions [6, 7] and beyond 
mean field studies with Gogny functionals [8, 9] have 
been performed. Mean field results in ^^S have shown 
the erosion of the = 28 shell closure and the manifes- 
tation of possible shape mixing and/or coexistence, ob- 
taining either rather flat potential energy surfaces (PES) 
or prolate-oblate minima along the axial quadrupole de- 
gree of freedom. Beyond mean field calculations -without 
particle number restoration- using either generator coor- 
dinate method combined with axial angular momentum 
projection or a collective hamiltonian have confirmed the 



relevance of allowing for configuration mixing in this nu- 
cleus. 

In this Rapid Communication we report the results 
of state-of-the-art Symmetry Conserving Configuration 
Mixing (SCCM) calculations to study the structure of 
^^S. The method is based on the Gogny DIS inter- 
action [10] and contains simultaneous particle number 
and angular momentum projection (PNAMP) of intrin- 
sic Hartree-Fock-Bogoliubov type states (HFB), includ- 
ing axial and triaxial shapes [11]. Variations of our ap- 
proach implemented with Skyrme [12] and Relativistic 
Mean Field [13] interactions have been also reported. 
The HFB-type states -|$(/32,7))- are found with the 
variation after particle number projection method (PN- 
VAP) [14, 15], being (/32,7) the quadrupole deformation 
parameters. Hence, contrary to the shell model, the col- 
lective intrinsic deformation is well established within 
this framework and this fact allows the description of 
the states in the laboratory frame in terms of their in- 
trinsic shapes unambiguously. The Ansatz for the wave 
functions is provided by the generator coordinate method 
(GCM) [15]. In this framework, the states are assumed 
to be linear combinations of particle number and angular 
momentum projected HFB-type states: 

l^/M.^^ E ^^^(/32,7)PifK^^i^^l^(/32,7)) (1) 

where /, M, K are the total angular momentum and its 
projection on the z-axis in the laboratory and intrinsic 
frame respectively, P^(^) the angular momen- 

tum and neutron (proton) projectors [15] and a labels dif- 
ferent states obtained for a given value of /. The param- 
eters gx{l^2il) are determined by the Ritz variational 
principle which leads to the Hill- Wheeler equation. The 
resulting matrices are diagonalized -one for each value of 
the angular momentum- providing the energy levels and 
collective wave functions defined in the (/32,7) plane. In 
addition, transition probabilities B{E2) and monopole 
strength p'^{EO) can be also calculated within this frame- 
work. Detailed expressions and performance of our ap- 
proach can be found in Ref. [11] (and references therein). 
This procedure upgrades in several ways the calculations 
performed in Refs. [8, 9]. On the one hand, intrinsic 
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FIG. 1. (color online) Potential energy surface -normalized to 
the absolute minimum- in the triaxial plane calculated with 
(a) PN-VAP approximation and (b)-(d) PNAMP approxima- 
tion with 7 = 0,2,4 respectively. Contour lines are separated 
by 1 MeV (dashed) and 2 MeV (sohd). 



HFB-type wave functions are found with the PN-VAP 
method instead of plain HFB improving the description 
of pairing correlations in the system [16, 17]. On the 
other hand, we include triaxial shapes -in Ref. [9] only 
axial configurations were allowed- and the GCM mixing 
is performed exactly without assuming a gaussian overlap 
approximation as in Ref. [8] . Finally, the absence of par- 
ticle number restoration leads to spurious mixing of solu- 
tions with different number of particles that could affect 
significantly the final spectrum [18]. Concerning techni- 
cal details about the calculations, we use a regular tri- 
angular mesh in the triaxial plane including A^gcm = 83 
HFB-type states. Each of these states is assumed to con- 
serve both time-reversal and spatial parity symmetries 
and is expanded in a single particle basis with nine ma- 
jor spherical harmonic oscillator shells. The number of 
integration points in the Euler -(a,6, c)- and gauge -ip- 
angles are chosen to ensure the convergence of both di- 
agonal and non-diagonal projected matrix elements of the 
total angular momentum and particle number operators 
to the nominal values /(/ + 1) and N = 28, Z = 16. In 
this case, these values are Na = 8, A/'^ = A/'c = 16 and 
N^ = 9 (see Ref. [11] for details). 

We first study the dependence of the energy of 
^^S with respect to the quadrupole deformation. In 
Fig. 1(a) the potential energy surface along the triax- 
ial plane calculated in the PN-VAP approximation is 
shown -the equivalent to the mean field energy in this 
method. We observe a very flat potential in a wide re- 
gion defined by the triangle which vertices are (/32,7) ~ 



{(0,0°), (0.4,0°), (0.3, 60°)}. For larger values of /^s we 
see a fast increase of the energy, having a slightly softer 
potential in the prolate part than in the oblate one. 
This behavior is a common feature of transitional nu- 
clei from spherical to deformed shapes. The minimum 
of the surface is not located at the spherical point but 
at (0.2,46°) which is an indication of the erosion of the 
N = 28 shell closure already at this stage of the cal- 
culations. To study the effect of angular momentum 
restoration in the system, we calculate at each point of 
the (/32,7) plane the particle number and angular mo- 
mentum projected (PNAMP) energy. In Fig. 1(b)- (d) 
we represent for / = 0,2,4 the corresponding PES, re- 
spectively. For / = we have only one projected PES 
while for / 7^ i^- mixing allows for up to //2 -|- 1 differ- 
ent solutions. Here, the lowest energy for a given value 
of the angular momentum at each (/32, 7) point is chosen. 
We now observe in Fig. 1(b)- (d) that the flat potential 
has been displaced to higher deformations excluding the 
configurations around the spherical shape that are now 
more than 4 MeV above the corresponding minima. Fur- 
thermore, we see that triaxial-oblate shapes are slightly 
more favored for / = around (32 ~ 0.3, while for / = 4 
triaxial-prolate shapes are a little more preferred around 
13 ~ 0.35. For the 7 = 2 case we observe an almost flat 
path in the whole range of 7 from (0.4,0°) to (0.3,60°). 
It is important to point out that we do not obtain PES 
with two (or more) well separated and almost degener- 
ated minima, which is one of the characteristic features of 
shape coexistence in nuclei -like, for example, in the pro- 
ton rich Se, Kr, Sr, Zr or Pb isotopes [19-21]. Strictly 
speaking we obtain flat potentials that suggest a wide 
configuration mixing rather than a shape coexistence. 

We now discuss the results of the SCCM calculations. 
The spectrum obtained after restoration of the symme- 
tries and shape mixing is plotted in Fig. 2(a). Here we 
have sorted the levels in bands, grouping states according 
to their largest B{E2) values. Nevertheless, we observe 
that there are relatively strong inter-band transitions 
showing that the underlying structure of these states has 
a pronounce mixing. The ground state is mainly popu- 
lated from the first 2^ excited state with an E2 transition 
and is also connected via a monopole transition {EO) to 
the O2" state. These two excited levels are very close in 
energy and they are also connected through an E2 tran- 
sition. The first 2^ state in the calculations is connected 
to several levels, mainly to the O2", 22 and 4^ but also 
to a lesser extent to 2^ and 43". This shows again the 
configuration mixing present in these states. 

To compare the results of our calculations to the 
available experimental information, we summarize in 
Fig. 2(b) the experimental level scheme of ^^S including 
all firmly established excited states Refs. [2, 4, 5]. 
The precise energy of the first excited 2+ state, 
E(2+)=1329(l) keV, has been determined in studies of 
the decay of the O2" isomeric state at 1365 keV [2, 4] in 
which the 2+ state decays at rest avoiding uncertainties 
due to the Doppler broadening of the lineshape. In 
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FIG. 2. (color online) Theoretical (left) and experimental [4, 5] (right) spectra for ^'^S. B[E2) values are in e^fm^, thick arrows 
represent ^0 transition with its corresponding value for . Experimental 7-rays and excitation energies are in keV. 



Ref. [5] three additional 7 rays with 949(5), 1128(6) and 
1929(7) keV have been observed in coincidence with 
the 2^ of ground state transition and assigned as 
directly populating the 2^ state from levels at 2278, 
2457, and 3248 keV, respectively. For the first two of 
these newly established states, spin values of 2+ and 4+ 
have been assigned on the basis of measured longitudinal 
momentum distributions [5] . A further 7-transition with 
an energy of 1891(10) keV, observed in coincidence with 
^^S residues, has tentatively been placed in the level 
scheme to connect the 3248 keV level to the isomeric 
O2" state based on the fact that the energy difference 
between the 1929(7) and 1891(10) keV 7-rays ( 38(12) 
keV ) is similar to the one between the O2" and 2^ 
states ( 36(2) keV ). A spin of 2+ has been tentatively 
assigned to the 3248 keV state due to its two decay 
branches to states with spin 2^ and O2". Finally, in 
the same work, a 2150 keV transition was tentatively 
placed as populating the ground state based on its 
non-observation in coincidence with the 1329 keV 2^ 
0]^ transition. However, this transition has been omitted 
in Fig. 2(b) since considering the quality of the spectra 
shown in Ref. [5] its placement seems much less justified 
as compared to the case of the 1891 keV transition, a 
judgement which is also supported by the huge difference 
between the experimental cross section for its population 
and the one calculated using reaction theory. The 
comparison between the experimental data and our 
theoretical results shows a nice qualitative agreement 
with a one-to-one correlation between observed and 



calculated states. In general the calculated excitation 
energies are larger than the experimental ones which 
can be explained by the lack of time-reversal symmetry 
breaking HFB states in the present work. Including 
those configurations we expect a compression of the 
theoretical spectrum [23]. Starting with the two first 
excited states at 1329 respectively 1365 keV, we obtain 
slightly larger values for both their energies as well as 
for the transition probabilities for their decay branches 
(see Fig. 2(b)). A similar picture is found in the shell 
model calculation of Ref. [4] -almost degenerated 2^ 
and O2" states and larger calculated B{E2) transition 
probabilities than experimentally obtained. However, 
the advantage of our approach as compared to the shell 
model is that we do not need to use effective charges 
and that we can also compute the value of the p'^{EO) 
within the same framework [22]. The latter quantity 
agrees remarkably well with experiment. 
Turning now to the higher-lying excited states, we will 
discuss in the following the main differences between 
the calculated and the experimental spectrum. For the 
22 state the calculations predict a second decay branch 
besides the observed 22 2^ transition, namely to 
the O2" state at 1365 keV. Such a transition, with an 
energy of 913 keV, has not been observed in Ref. [5] 
and this non-observation is significant given the high 
statistics obtained in that experiment. The calculated 
4^ state decays mainly to the 2f level as is observed 
experimentally. In addition, a second decay branch to 
the 22 state is obtained which has not been observed 
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experimentally. However, since experimentally the 4^ 
state is observed at much lower excitation energy, this 
low-energy 7-transition of 179 keV is hindered due 
to phase-space arguments and furthermore may have 
escaped observation due to the relatively high energy 
thresholds used in the experiments. 
Also in the case of the third excited 2+ state, the 
calculations predict one more decay branch than iden- 
tified experimentally. While the calculated decays of 
the 23" state to the 2^ and O2" levels can be identified 
with the 1929(7) and 1891(10) keV 7-rays observed in 
Ref. [5], a 980 keV 2^ 22 transition has not been 
reported. But again, the transition is hindered due to 
the lower transition energy as compared to the other two 
decay branches which may explain its non-observation. 
According to our results, the 2^ state is the band head 
of a quasi-7 band with levels 3^,44 on top it. 

To shed more light on the structure of the calculated 
states we represent in Fig. 3 the collective wave functions 
of the most relevant ones. Additionally, their distribu- 
tion in K and spectroscopic quadrupole moments are also 
given in Table I. In Fig. 3 we observe that the distribution 
of probability in the (/32,7) plane for each state belong- 
ing to the same band changes significantly. However, 
shape mixing allows the overlap between these states 
and the quadrupole transitions between them. Hence, 
the ground state 0^ is distributed almost uniformly in 
the whole range of 7 in a trapezoidal region between 
(0.1, 0°), (0.3, 0°), (0.1, 60°), (0.35, 60°). Therefore, the 
ground state for ^^S does not correspond to a spheri- 
cal shape -as it should be if the N = 28 shell closure is 
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FIG. 3. (color online) Collective wave functions in the triaxial 
plane for the states represented in Fig. 2. 



TABLE I. Distribution of < 4 and spectroscopic 

quadrupole moments of the states appearing in Figure 3. 



preserved- but a mixing of deformations where the triax- 
ial degree of freedom plays an important role. The rest of 
the states of this band, which have a = configuration, 
evolve towards axial prolate shapes. Their maxima are 
found at /32 ^ 0.35, having the 2^ more mixing with tri- 
axial and axial oblate configurations than the 4^. These 
prolate states appear in shell model calculations [4] al- 
though in that case a prolate shape is assumed for the 
Oi contrary to the result found here. Furthermore, in 
Ref. [5] the 4f state is predicted to have a neutron Ip-lh 
excitation configuration, being the band head of a K = 4 
band. This is not the case in the present work where the 
first 4+ state with the highest K = 4 component is found 
at 5.4 MeV having a small B{E2) value to the 2f state. 
The first excited band is also K = mainly. The band 
head O2" state is axial prolate deformed with the max- 
imum of the probability found at /32 ^ 0.4. This dis- 
tribution explains the large transition probability to the 
2^ state. We also observe an axial prolate-oblate mixing 
for the 22 state which connects to the prolate deformed 
O2", 2^ and 4^ states and also to the mainly oblate con- 
figuration of the 42 state. A rather small spectroscopic 
quadrupole moment of the 22 state is found here and 
also in shell model calculations. This result was inter- 
preted as an indication of spherical shapes for both 22 
and O2" [4]. However, this is not the case in the present 
calculations. 

Finally, a quasi-gamma band with K = 0^2 mixing for 
/-even and i^' = 2 for /-odd is obtained grouping the 2^ , 
3i and 4^ states. The collective wave functions of these 
states are rather similar, with the corresponding max- 
ima found at ^ (0.3,40°) and the probability decreasing 
slightly more pronouncedly in the 7 direction than in the 
P2 one. The values for the spectroscopic quadrupole mo- 
ments are typical of intrinsic oblate shapes when K = 2 
configurations are also present, i.e., negative, zero and 
positive for the 2+,3+,4+ states respectively. 
To summarize, we have study the nuclear structure of ^^S 
with a state-of-the-art SCCM method including beyond- 
mean-field effects such as particle number and angular 
momentum restoration and shape mixing of axial and 
triaxial configurations. Wide configuration mixing rather 
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than shape coexistence is found in most of the low-lying 
states supporting the erosion of the TV = 28 shell clo- 
sure and the transitional character of this isotope. A 
deformed and 7-soft 0^ as well as a prolate deformed O2" 
states are obtained. Our results are in a very good quali- 
tative agreement with the recently available experimental 
data and with the most recent shell model calculations 
although the interpretation of some of these data is dif- 
ferent. Finally, we also find a quasi-7 band with decaying 
properties that could be attributed to the experimental 
(2+) level measured at 3248 keV. 
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